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Abstract. We consider the relation between exact solutions of cosmological 
models having minimally and non-minimally coupled scalar fields. This is done for 
a particular class of solvable models which, in the Einstein frame, have potentials 
depending on hyperbolic functions and in the Jordan frame, where the non- 
minimal coupling is conformal, possess a relatively simple dynamics. We show 
that a particular model in this class can be generalized to the cases of closed 
and open Friedmann universes and still exhibits a simple dynamics. Further we 
illustrate the conditions for the existences of bounces in some sub-classes of the 
set of integrable models we have considered. 
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1. Introduction 

Recent observations [HE] support interest in inflationary models with non-minimal 
coupling between the inflaton field and the Ricci scalar curvature R. There are two 
classes of non-minimally coupled inflationary models. In one of them a function of 
the scalar field (without a constant term) non-minimally coupled to R is added to the 
Hilbert-Einstein term while in the other it replaces the latter (4). This second 
class of models is known under the name of induced gravity and cosmology. 

In spite of great success of numerical simulations and of using of different 
approximation schemes for studying of cosmological solutions of the Einstein 
equations, the exact analytical solutions of these equations are always useful. Indeed, 
they permit to catch some important qualitative features of different cosmological 
models in the most easy and convenient way. 

Most of the results of the exact integration of cosmological models with scalar 
fields are connected with minimally coupled models. A rather exhaustive list of such 
models was presented in paper [^. 

The integrability of the models with minimal and non-minimal coupling was 
studied by using some advanced mathematical methods in [^. 

In |7] we proposed a method for constructing integrable models with non- 
minimally coupled scalar fields by using the interrelation between the Jordan and 
Einstein frames. We considered some explicitly integrable cosmological models 
for spatially flat Friedmann universes filled with a minimally coupled scalar field, 
introduced in (H. On using the Einstein-Jordan frame transition, we constructed 
their counterparts in induced gravity models and in models having a conformally 
coupled scalar field plus the Hilbert-Einstein term. 

In a recent paper |8] , a flat Friedmann cosmological model in which the scalar field 
was conformally coupled to gravity and its potential included a positive cosmological 
term and a negative term with a quartic self-interaction was studied. Such a model 
is exactly integrable and for a large class of initial conditions possesses a bounce 
and avoids the cosmological singularity. Further, in such a model the cosmological 
evolution is, to a large extent, independent of the evolution of the scalar field. Let us 
note, that the presence of bounces was considered as a positive feature in cosmology 
and their study has rather a rich history [^ . In a subsequent paper [lOj , the minimally 
coupled counterpart of the above model was considered. The corresponding potential 
is a linear combination of hyperbolic sine and cosine functions to the fourth power. It 
is curious that the model with such a potential was previously studied in great detail 
in a series of papers [nKHEHIKIi]. In these papers such interesting topics as the 
possibility of the crossing of the Big Bang - Big Crunch singularities and the existence 
of anti-gravity regimes were studied. On the other hand one can see that this minimally 
coupled model belongs to a more wide class of models, described in the paper |5| as 
a class 7. We did not consider the non-minimally coupled counterpart of the models 
of this interesting class in our preceding paper |7] and in the present paper we want 
to fill this gap. Hence, we study the models of this class and find the conditions 
whereby they possesses a bounce. Further, we show that the model, considered in 
mm is easily generalizable to the cases of closed and open Friedmann universes, and 
present the corresponding solutions. Remarkably, in this case the integrability of the 
non-minimally coupled model is more apparent that the integrability of its minimally 
coupled counterpart. 

The structure of the paper is as follows: in Sect. 2 we review formulae that 
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connect models in the Einstein and Jordan frames; in Sect. 3 we describe the integrable 
cosmological model with a conformally coupled scalar field and its generalisation for 
the cases of closed and open Friedmann universes; in Sect. 4 we discuss the properties 
of the cosmological solutions for a general class of minimally coupled models with 
potentials, depending on hyperbolic functions and for their non-minimally coupled 
counterparts; in Sect. 5 we study when the bounces are possible in these conformally 
coupled models; in Sect. 6 we write down some formulae the corresponding induced 
gravity models. The last Section contains concluding remarks. 


2. Relations between general exact solutions in the models with 
minimally coupled and non-minimally coupled scalar fields 


Let us consider a cosmological model, described by the following action 


S= d‘^Xy/—g 


U{(t)R - - V(a) 


( 1 ) 


where U{a) and V(a) are differentiable functions of the scalar field a. If U{a) is a 
constant: U{a) = Uq, then the scalar field is minimally coupled to gravity, otherwise 
we have a model with non-minimal coupling. 

Let us consider the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric. In 
this case the interval is 


ds^ = NirYdr^ — a^Y 


l-Kr^ 


— r^dO^ — s\rr(9)dip‘ 


( 2 ) 


where a(T) is the scale factor and N(t) is the lapse function, K is a constant. As 
usual K = 0 describes a flat universe, AT = 1 a closed universe, and AT = —1 an open 
one. On substituting this metric into the action o and varying it with respect to IV, 
a and cr, we obtain two Friedmann equations and a Klein-Gordon equation: 


6t7 


,2 ' 
h H- 9 “^ 


+ QhU'a = + N^V, 


2U 







2h + 3h^ - 2h— + —K 

+ 2U' 

a + 2ha - a 

= N'^V - 

2U" + 1 

N a? 

N 


2 


( 3 ) 

( 4 ) 


( 


n N 

3d — — 

a-6U' 

h + 2h^ -h— + —K 

V 


N 


+ iVV' = 0, (5) 


where h = a/a, a “dot” means a derivative with respect to the parametric time r and 
a “prime” means a derivative with respect to the scalar field. 

If U{a) is not a constant, then the scalar field a is non-minimally coupled to 
gravity. Let us perform a conformal transformation of the metric 
Uo 


where Uo is a positive constant. We also introduce a new scalar field (p such that 


( 6 ) 


_ y'UojU + iU'^) 

da U 




U 


da. 


( 7 ) 


As a result, the action © becomes the action for a minimally coupled scalar field: 


S = 


d‘^X'Y—9 


UoRi'g) - 


( 8 ) 
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where 


W{^) = 




Let us emphasise that the formulae (jS])-® are valid for any metric. 
In the Einstein frame the FLRW metric (HI) becomes 

dr^ 


ds^ = - N^dr^ + 


— s,iir{9)dip" 


, 1 - Kr"^ 

where the new lapse function and the new scale factor are defined as 



For the case K = Q the Friedmann equations in the Einstein frame are: 

4Uoh + eUo'h'^ - iUoh^ = + N^W, 

N ^ 


0 + 



^ + N^w;^ = o, 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


where ft, = d/d. 

Our goal is to obtain the general solutions for the Friedmann equations both in 
the Einstein and Jordan frames. We are assuming that such solutions are known in 
one frame and shall seek them in the other one. The algorithm has been proposed in 
our paper [7]. In this paper we consider new interesting examples and generalize our 
algorithm to the cases of open and closed Friedmann universes. 

In [7] we considered in detail two specific forms of the function U{a). For the 
induced gravity case mile] 


Uind{^) — 2 ^ 


then we have from Eq. dZD 


/2t7o(l + 67 ) 


In 


7 


a 

CTo 


and, inversely, cr = 


We assume that 7 ^ —1/6 since for the case of conformal coupling, 7 
nontrivial solutions can exist only for the potential V = m- For 7 = 
add a positive constant Uq to the induced gravity term and consider [7]: 


Uc{a) = Uo-j^, 


(15) 


(16) 

= - 1 / 6 , 

— 1/6 we 


(17) 


i.e. the case for which the coupling is conformal and a nonzero Einstein-Hilbert term 
is also present [suniiin]. For this case we get 


(j) = \/3t/o In 


yml^+a 

vOTo-d 


and 


d = \/l2Uo tanh 




(18) 
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The model with U{a) = Uc and the potential 

K = A - ccr^ , (19) 

with A > 0 and c > 0 has been considered in [8]. It was shown that the first 
Friedmann equation for this model reduces to the standard Friedmann equation 
where matter is represented by two perfect fluids: a positive cosmological constant 
and radiation, whose energy density can be both, positive or negative, depending 
on the initial conditions for the scalar field. If the energy density of the effective 
radiation is negative, then the universe undergoes a bounce. In the subsequent paper 
m, it was shown that on applying a conformal transformation, combined with a 
reparametrisation of the scalar field (just as the procedure described in our paper [7] 
and in the preceding section of the present paper), one arrives to a minimally coupled 
model with a potential 


W{(j)) = A cosh"*^ 


- IddC/p^csinh^ 


( 20 ) 


2y/w^j " V2V3C7o/ 

As we have already mentioned in the Introduction, this model was intensively 
investigated in papers [IllllllllsllIlllIS]. In particular, the corresponding inflationary 
scenario has been considered in m- 

Further, it belongs to the set of integrable models [5] with a minimal coupling 
and a potential 


W{(f)) = Cl I cosh 


3/3 


2(l-/3) 




, _ , 1 C 2 I sinh ,_ 

.ymToj; V LyisiTo. 

where ci, and C 2 are arbitrary constants. The model with the potential dini) 
corresponds to the particular case of the potential m with /3 = 1/3. We shall discuss 
some properties of the above family of models in the subsequent sections. Here we 
wish to stress that the main attractive feature of this class of models is their exact 
integrability in terms of the cosmic time parameter. Thus, it gives the opportunity of 
studying such fundamental questions of cosmology as the behavior of the universe in 
the vicinity of the cosmological singularity and the relations between different frames 
and different parameterizations of scalar fields. The question whether the integrable 
one-field cosmological model classified in can be embedded as consistent one-field 
truncations into Extended Gauged Supergravity or in A = 1 supergravity gauged by 
a superpotential has been considered in |20| . 

We now present the generalisation of the model of paper [8] to the case of closed 
and open Friedmann universes. 

The first Friedmann equation ([3]) for this model is 


3/3(/ 


2 ( 1 - 3 ) 

3 


,( 21 ) 


6 t/n- 


+ 6 Go - 


1 

12 ^ 

The Klein-Gordon equation ([5|) is 


1 

12 ^ 


KN^ 


1 


+ N‘^{A-ca^). 



a -\- a 


a a“ \ dN , o i 

2 I aN 


a" 


N^Ka 


= 0 . 


( 22 ) 


(23) 


On making a substitution such as the one done in paper [5]: cr = y/a, in the 
Klein-Gordon equation (l23l) . we reduce it to 

xN dcy^ ^ N'^Kx 


.. , Xa 
x +- 


= 0 . 


(24) 
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On choosing the lapse function N as N = a, we can rewrite Eq. (IMl) in terms of the 
derivatives with respect to the conformal time rj: 


^ - 4cx^ + Kx = 0- 

arj^ 


(25) 


On multiplying this equation by the conformal time derivative of x, we find its 
first integral 

2 


^ f ^ ) -cx^ + = A 


2 \ dr] 


1 


(26) 


where >1 is a constant. Let us now rewrite the expression p 6 p in terms of the initial 
scalar field a and time derivatives in terms of the cosmic time (i.e. to the time 
parameter corresponding to the lapse function choice TV = 1). One obtains 

ifcr2 A 


1 • 2 I 1 2 • 2 

-CT + -a a 


, crcr- - ccr + 2 

2 2 a 2a‘‘ 

On comparing Eq. (I27p with the Friedmann equation (1^^ . we see that the latter 
reduces to a very simple form 

6Uq—+ 6Uo—= A-\ —j. (28) 

This is a natural generalisation of the Friedmann equation, obtained in paper [ 8 ] . The 
equation for possible turning points (bounces or the points of maximal expansion) is 


(27) 


a^- 


6UoK 2 A 


A 


+ -- 0 . 


(29) 


Let us give a list of the cosmological evolutions for different choices of the curvature 
K and the radiation constant A. 

1. (a). K = 0, A < 0. We have a bounce at 


as = 



This case was analysed in detail in papers 

(b) K = 0, A = 0. For this case we have an infinite de Sitter expansion or an infinite 
de Sitter contraction. The singularity is absent. 

(c) K = 0, A > 0. We have an infinite expansion which begins from the Big Bang 
singularity or an infinite contraction which ends in the Big Crunch singularity. 

2. (a) K = 1, A < 0. We have a bounce at 


Ob = 




(h) K = 1, A = 0. We have a closed de Sitter universe which contracts, has a bounce 
at 


Ub = 
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and then expands infinitely. 

(c) if = 1, 0 < For this case we have two types of possible evolutions. One 

begins at the Big Bang singularity, expands until the point of maximal expansion 


Om 




and then contracts until the encounter with the Big Crunch singularity. The second 
type of evolution is that with the bounce at 



(d) K = 1, A = We have an Einstein static universe with radius 

3t/o 

a5 = —. 

9C/^ 

(e) K = 1, A > —We have an infinitely expanding or an infinitely contracting 
universe. 

3. (a) K = —1, A < 0. We have a bounce at 


Ob = 


A V A2 a) 


(b) K = —1, A ^ 0. We have an infinite expansion beginning from the Big Bang or 
an infinite contraction culminating in the Big Crunch. 

For all these cases the Friedmann equation (l2^ can be integrated explicitly, because 
the Ricci scalar as an integral of motion. Combining Eqs. ([S])-® it is easy to show 
that R = Rq = 2A/Uq at all values of K. We shall write down the explicit solution 
only for the case 2 (a). The expression for the scale factor is 

1/2 


ait) = 




A 9 “ T 

A2 A 


A 


24t/f 


■t 


(30) 


The Hubble parameter is 


H{t) = 


A\/9C/^ - AAsinh 


96Uq 


3C/o + v'9C/2-AAcosh 


(31) 


Remarkably all the dependence of the cosmological evolution on the scalar field 
is encoded in the quantity A, i.e. it is nothing more than the evolution of the universe 
filled with the cosmological constant and the radiation-type fluid (see Eq. (I^ b 

Let us notice that this “disentanglement” of the cosmological dynamics from the 
dynamics of the scalar field is a rather unusual phenomenon. For the case of models 
with a minimally coupled scalar field the only potential where this takes place is the 
constant potential, i.e. a cosmological constant. For the case of a conformally coupled 
scalar field the simple form of this Friedmann equation, which reduces to the standard 
one for a mixture of fluids, is connected with the fact that the conformally coupled 
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scalar field behaves just like a radiation fluid with the constant A, which exhibits 
the properties of such a fluid. It is also important that quartic self-interaction term 
is also conformally invariant. The presence of this term gives us the opportunity of 
making the constant A negative, by choosing a negative sign for the constant of the 
self-interaction of the scalar field. 


4. General solutions for a class of minimally coupled models with 
potentials, depending on hyperbolic functions and for their 
non-minimally coupled counterparts 

The integrability of the class of models with the potential (l?T|) was discussed in 
We shall arrive to this result in a slightly different way and also write down the 
corresponding formulae for the conformally coupled counterpart of this class. 

The minisuperspace Lagrangian, generating the Friedmann equations ca), (USD 
and the Klein-Gordon equation m is 

6(a)^a[/o 


L = 


NWd^. 


N 2N 

If one considers the potential (|2T]) and chooses the lapse functions 
N = 1^53-6/3 

then the Lagrangian has the following form: 

9 417o 


(32) 


L = L/ 32 --- 6 / 3-2 _ 
2 


8C/0 




Let us introduce new variables x and y defined as 

/ 6/3 A X 


= xy, 


exp 


then the Lagrangian ([M)l takes the form 




-w. 


(33) 

(34) 

(35) 


For the potential W, given by m, we finally obtain, on introducing another 
couple of independent variables 

X + y X — y 




V = 


2 2 
the following simple expression for the Lagrangian: 

. 417o 

; 3/32 


L = 


(ciC 


2(l-/3) 


2 (l-g) 


C 2 ?? 


The corresponding Euler-Lagrange equations are 


.. 417o^2(l-^) 


3/32 


2 Q 

e?-3 = o. 


4C/o^2(l-/3) 


Their first integrals are 


^ 4C/o 1(1^ _ 

2 ^ 3 / 32 ^ 


al 

2 


C 2 


/3 

417o 


2 o 

^ = 0 . 


(36) 


(37) 


(38) 


3/3- 




= E2. 


(39) 
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These equations have solutions by quadrature. For /3 = 1/3 one obtains the 
solutions in terms of elliptic functions m- 

The case /3 = 1/2 is also well-known. For this case the parametric time coincides 
with the cosmic time and the solutions can be expressed in terms of elementary 
functions (see e.g. [21] and references therein). 

It immediately follows from the first Friedmann equation that Ei = i?2. We can 
now find the cosmological variables 


= - n 


2\1/(6/3) 


N = 




(l-2/3)/(2/3) 




One can also write down the expressions for the solutions of the corresponding 
models with non-minimal coupling. Indeed, on using (HU cLnd (j ISp . w6 obtain 


«. = p (({++ («-1)'/'’"’) {e - , 




for the model with conformal coupling. At ,5 = 1/3 we get Qc = 
The associated explicit expressions for the potential 14 is 

1 f [(yT2C7o + a)3/3 + (vOTo-a)3/5]^ ^ 
3617241//3 I (I2C/0 - a2)i-3/3 

[(v^T2C7o + - (vOTo - a)3/5] ^ \ 

^ """ (12C/o - u2)l-3/3 f ■ 


(40) 

(41) 

(42) 


(43) 


and corresponds to the potential m- In the next section we shall consider some 
general properties of such potentials. 


5. Integrable models with bounce solutions 

Let us consider the model with a conformal coupling, the Hilbert-Einstein term and 
the potential 14, given by (P|. We wish to find the conditions for the existence 
of bounces in such models. Let us first consider the Friedmann equations and the 
Klein-Gordon equation for an arbitrary model, wherein the scalar field is conformally 
coupled to the scalar curvature. 

A bounce point ts is defined by two conditions: the Hubble parameter H = 0 
and H > 0. On substituting = 1 and H = h = 0 into system (jSj)-® with K = 0, 
we get 

l^^+V = 0, (44) 

AUH + ( 2 U"+ + 2U'd-V = 0, (45) 


a - QU'H + V' = Q. 


(46) 
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On substituting U = Uc and using (Hill and (H51) . we obtain from Eq. (H5|) 

4UoH-^V+^aV' = 0. (47) 

From Eq. (SZl) one has 

i7=^(4E-uE0. (48) 

Thus, we have the condition 

4V - aV' > 0. (49) 

Further, it follows from Eq. (|4^ that at the point where H — 0, the potential E < 0. 

We can now consider some particular values of the parameter /S, for which the 
potential (|^ appears relatively simple. 

For the case (3 = 1/3, the potential (l4^ can be written as 

E = El + Eacj^ (50) 

On substituting formula (l50l) into formula (|4^ we obtain the condition 

El > 0. 

Further, in order to have a negative potential (150)1 at some value of the scalar field it 
is necessary to have E2 < 0. This is just the case considered in [Hilo]. 

For the case /? = 1, the potential (1431) can be written as 

V = V^{uo-'^ , (51) 

where Vq = (ci + C2)/17o- This potential can be negative if and only if the constant 
Eo < 0. Then, on substituting the potential (IQTl) into the condition (HOI) , we obtain 

12C/o < cr^ (52) 

Another interesting value is /3 = 2/3. In this case the potential can be written 
as 

= Y^(12C/o-a2)(a2 + 12[/o + 2Eia), (53) 

where Ei = 2c2-\/3C7o/ci. The study of the conditions for the existence of bounces 
with arbitrary values of the constants ci and Ei is rather cumbersome. Thus, we shall 
limit ourselves to considering a few particular values of the constant Ei. 

Firstly, let us consider the case, when 

El = - 

Now, the potential E can be written down as 

while 

4E - ctE' = (a - y^l^) " (a + 2y^l^) . 

In this case two conditions of the existence of bounce are compatible if 
Cl > 0, a> v^l^ or — 2y'l2Uo <a< -Vl2[/o. 
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Secondly, let us consider the case 
Vi = v^T^. 

Now the potential is 


V = 


Cl 


U4[/g 


( 12 Uo-ct^) (^a+^/T^y 


and 


4y-at/'= -2\/l^^^^(a-2v^l^) (fT+V^l^) . 


The conditions for the existence of bounces are compatible if 

Cl > 0, v'121/o < tr < 2i/l2Uo. (54) 

We see that in these cases as well as for /3 = 1 the bounce corresponds to a such 
tr that U < 0. The following example shows that for U > 0 the bounce is possible not 
only in the case /3 = 1/3. Let 


289^2 

For positive values of a we get that F < 0 at 
— i/QUo < a < i/lWo- 


Vi = -^-^Vm. 


(55) 


(56) 


At Uq = 1, j^-\/6t/o ~ 0.20412. At the same time the condition (l4^ is satisfied if 
a < 0.27178 (see Fig. [T]). Thus, we get a bounce point at 17 > 0. 



Figure 1. The potential Vc (solid red line) and the expression AVc —crVy (dash 
blue line) at Uq = c\ = p = 2/3, Vi = —289\/6/24. The fine structure, 
presented on the right picture, shows that a bounce is possible. 


For the case /3 = 4/3 the potential can be represented as: 

(l2{7o - cr^)^ C 2 4ci 

576t7o' \/2 ^4 +144 [/oa2 + 288170/ 

In this case it is difficult to obtain solutions in terms of elementary functions. We 
consider it by using numerical calculations. For 17o = 1/4, ci = —1, and C2 = 0.7 we 
demonstrate (see Fig. [2]) that a bounce is possible at U{a) >0 and the corresponding 
potential has a minimum. 

Numerical calculations show that in this case the conditions 14 < 0 and 
4Vc - crV/ > 0, require cr G (0.230742, 0.285248). 
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Figure 2. The potential Vc (solid red line) and the expression AVc —crV^ (dash 
blue line) and the function U (dot black line) at Uq = 1/4, ci = —1, C 2 = 0.7, and 
0 = 4/3. The fine structure, presented on the right picture, shows that a bounce 
is possible. 


6. The induced gravity models 


On using ([HI) and (fT6|l . we construct the induced gravity model 
S' = y {UindR - + Vind'^ , 


(57) 


with the following potential 


^ind — 




4^ C/2 




Cl I + cr. 




2(l-/3) 


+ C 2 cr 


(‘ 


.6/3r_^6/3r)- 


2(l-/3) 


where T = 

In order to obtain the general solution we use formulae (HU and (HU and write 
the general solution in terms of the variables ^ and rj: 

i/(6/3r) 


^ind — ^0 


i±^\ 


O^ind — 


O-Q 

N.nd = 

3p2cro 


C + r?y 


i/(6/3r) 


ie-v^) 


l/(6/3) 


l/(6/3r) 




) 


ie - -t) 


-l + l/(2/3) 


(58) 

(59) 

(60) 


Thus, in order to obtain the general solutions for the induced gravity model (l57l) . 
one has to solve the same equations (l38l) as for the model with a minimally coupled 
scalar field and the potential m- Note that for /3 = 1 /3 this induced gravity model 
is integrable for nonzero values of K as well. 


7. Concluding remark 

On using the connections between exact solutions of the cosmological models with 
minimally coupled and non-minimally coupled scalar fields, one can obtain some 
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interesting results. We have applied this connection to show that the integrable 
cosmological model [5] with a constant R belongs to a one-parameter set of integrable 
models with the same function Uc and different potentials 14. The corresponding 
integrable induced gravity models were also constructed. 

A cosmological model, with a particular potential depending on hyperbolic 
functions, was intensively studied in a series of papers [in na na [H [15] having 
a special emphasis on questions such as the possibility of crossing the Big Bang- 
Big Crunch singularity and an antigravity regime. Such considerations led to some 
discussions Ha HSl IMl Hi]- Remarkably, this model has a conformally coupled 
counterpart in an integrable model with a rather simple cosmological dynamics and a 
bounce solution, thus avoiding the cosmological singularity laiio]- 

On the other hand, this model with its potential, which includes fourth degrees 
of hyperbolic functions, can be considered as a particular case of a wider class of 
models (case 7 in paper 0), for which, in the present paper, we have constructed 
their conformally coupled and induced gravity analogues. We note that sometimes 
the model with non-minimal coupling can be integrated on using some symmetry 
properties (see, for example [26] 1. so it may be useful to consider conformally 
connected models both with minimal and non-minimal couplings in order to obtain 
new integrable models. In this paper we find integrable cosmological models with 
non-minimally coupled scalar fields with both open and closed FLRW metrics. 

Let us emphasize that the models with a minimal coupling cannot have bouncing 
solutions since the second time derivative of the scale factor is proportional to the 
time derivative squared of the scalar field taken with the opposite sign. Thus, the 
dynamics in models with a minimal coupling and that in their non-minimally coupled 
counterparts can be very different. We shall study this question in detail in a future 
publication Ezj. Finally, we wish to end by noting that the relation between the exact 
solutions in non-minimally coupled and minimally coupled cosmological models can 
be used to address some additional points in discussions concerning the crossing of 
the singularity m- 
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